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ABSTRACT 


An  analysis  of  a  doubly  coupled  single  fiber 
resonator  has  been  made.  The  addition  of  a  second 
coupler  to  a  typical  single  coupler  resonator  with  less 
than  100  percent  coupling  is  shown  to  greatly  improve  the 
finesse  and  sensitivity  to  fiber  attenuation  of  the 
resonator.  A  simple  method  of  achieving  extremely  high 
finesse  with  approximately  70  percent  coupling  is  also 
shown.  With  typical  fiber  and  coupler  loss  values  from 
the  literature  included  in  the  calculations,  this  device 
shows  a  minimum  detectable  signal  on  the  order  of  20 
picometers  as  a  microbend  sensor  and  .17  microradians  per 
second  as  a  laser  gyroscope. __ 
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ABSTRACT 


An  analysis  o£  a  doubly  coupled  single  fiber 
resonator  has  been  made.  The  addition  of  a  second 
coupler  to  a  typical  single  coupler  resonator  with  less 
than  100  percent  coupling  is  shown  to  greatly  improve  the 
finesse  and  sensitivity  to  fiber  attenuation  of  the 
resonator.  A  simple  method  of  achieving  extremely  high 
finesse  with  approximately  70  percent  coupling  is  also 
shown.  With  typical  fiber  and  coupler  loss  values  from 
the  literature  included  in  the  calculations,  this  device 
shows  a  minimum  detectable  signal  on  the  order  of  20 
picometers  as  a  microbend  sensor  and  .17  microradians  per 
second  as  a  laser  gyroscope. 
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I.  INTRODUCTION 

The  use  of  evanescent  wave  coupling  in  opti¬ 
cal  fibers  allows  for  the  design  and  implementation  of 
elaborate  feedback  networks.  Since  reflections  do  not 
occur  at  the  coupling  interface,  the  analysis  of  these 
devices  is  quite  straightforward  and  leads  to  some  very 
useful  results. 

One  of  the  dominant  applications  of  fiber 
optic  technology  is  in  the  sensor  research  area.  Most 
fiber  optic  sensors  can  be  divided  into  two  broad  catego¬ 
ries:  those  that  employ  an  amplitude  modulation  mecha¬ 
nism  for  detection  and  those  that  employ  a  phase  modula¬ 
tion  mechanism  [3].  Application  of  these  mechanisms  to 
sensor  devices  include  interferometric  acoustic  sensors 
[2],  fiber  microbend  acoustic  sensors  [18],  fiber  optic 
accelerometers  [19],  magnetic  field  sensors  [20],  and 
fiber  optic  gyroscopes  [14,  15].  Another  rapidly  growing 
application  area  is  signal  processing  [12],  specifically 
the  use  of  fibers  to  form  a  recirculating  delay  line  [10] 
and  the  use  of  an  evanescently  coupled  feedback  network 
to  perform  matrix  multiplication  [11]. 
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The  purpose  of  this  thesis  is  to  analyze  a 
particular  resonant  feedback  structure  and  describe  its 
unique  characteristics.  Application  of  this  device  as  a 
microbend  sensor  will  be  discussed  in  detail  and  sugges¬ 
tions  made  for  other  applications. 

II.  SIMILAR  STRUCTURES 

In  general,  resonators  are  extremely  sensi¬ 
tive  to  losses  in  the  feedback  path.  Since  low  loss 
fiber  is  readily  available,  the  key  component  in  the 
development  of  a  high  finesse  fiber  optic  resonator  is  a 
low  loss  coupler.  In  reference  16],  a  mechanical  lapping 
technique  for  the  fabrication  of  a  single-mode  fiber 
optic  directional  coupler  was  introduced.  The  coupling 
in  this  device  was  tunable  and  achieved  a  power  transfer 
of  70  percent  with  less  than  10  percent  loss.  When  this 
coupler  was  first  used  in  an  all-single-mode  fiber  reso¬ 
nator  18],  a  finesse  of  80  was  achieved.  The  dominant 
loss  in  this  structure  was  the  coupler  insertion  loss, 
approximately  3.2  percent. 

In  ref.  (7],  a  refined  mechanical  lapping 
technique  was  used  to  fabricate  a  low  loss  coupler.  This 
tunable  device  coupled  up  to  99.8  percent  of  the  input 
power  with  losses  approaching  .1  percent.  In  addition. 
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the  authors  found  that  the  scattering  loss  varied  with 
the  coupled  power.  This  low  loss  coupler  was  used  to 
construct  a  low  loss  fiber  ring  resonator  [9].  For  a  one 
meter  loop  loss  of  .41  percent,  the  authors  measured  a 
finesse  of  approximately  500  with  a  coupler  loss  of  .22 
percent.  Therefore,  for  the  first  time,  the  resonator 
loss  was  dominated  by  the  fiber  loss  rather  than  the 
coupler  loss. 

The  development  of  these  low  loss  fiber  reso¬ 
nators  will  open  the  door  to  a  whole  new  family  of  sen¬ 
sors  and  optical  processing  devices.  These  published 
results  will  be  applied  to  the  derivations  in  the  appli¬ 
cations  section.  This  will  allow  for  a  realistic  simula¬ 
tion  of  the  device's  operation. 

III.  GENERAL  DEVICE  DESCRIPTION 

The  main  feedback  structure  of  this  device 
can  be  constructed  from  a  single  length  of  optical  fiber. 
A  general  diagram  is  shown  in  Figure  1.  The  supplemental 
paths  A  and  5  can  be  implemented  by  using  separate 
lengths  of  fiber  or  placing  a  loop  in  the  main  fiber. 
Note  that  this  design  reduces  to  a  simple  resonator, 
described  in  Section  II,  when  there  is  zerc  coupling  at 
CPI  and  the  supplemental  paths  are  removed. 
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Using  coupled  mode  theory  [1]#  the  evanescent 
wave  coupling  that  occurs  at  couplers  one  and  two,  CPI 
and  CP2  respectively,  can  be  described  by  the  following 
equations : 


E3  * 

■  ^^1^1  ^1®2 

E?  »  t2E5  +  ir2Eg  (3) 

Eg  »  ^^^2^5  ^2®6 

r^  +  t^  +  Y*  1  (5) 

where:  -  coupled  power 

t^  ®  through  put  power 

Y  =  coupler  power  loss 


The  feedback  loop  and  supplemental  paths  provide  the 
following  additional  relationships: 


E5  ® 

E3 

(6) 

^6  = 

(7) 

E2  ® 

X3e^PE-7 

(8) 

where: 

«1 

X2 

“2 

X3 

P 


transmission  coefficient  of  path  A 
phase  shift  in  path  A 
transmission  coefficient  of  path  B 
phase  shift  in  path  B 

transmission  coefficient  of  feedback  loop 
phase  shift  in  feedback  loop 


■  -.'‘K 


‘•■V*  ' 


M 
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Using  equations  1-8,  we  can  determine  the  value  at 

any 

port  with  respect  to  the  input.  These  values  are  deter- 

mined  by 

the  following  equations: 

NOTE; 

these  equations  assume  that  r  and  t  satisfy 

equation  5. 

d  =  1  - 

XiX3rit2e^<P+“i+''/2)  _  x2X3r2tiei<P+°‘2+’^/2) 

(9) 

Eo  =  1 

^  a 

(XiX3tit2e^^P'‘'“i^  -  X2X3rir2e^^P'*‘“2^  )E]^ 

(10) 

1 

i(p+a2+iT/2) 

^3  =  d 

(ti  -  X2X3rir2e  )  E^ 

(11) 

^4  =  Z 

(r^e^^/^  +  XiX3t2e^^P‘'‘“i^  ) 

(12) 

^5  =  Z 

(Xitj^e^^i  -  X3^X2X3r2e^^P'^“i'^“2'*’'^/^^  )  E^ 

(13) 

^6  =  Z 

(X2rie^^°‘2+Tr/2)  +  Xj^X2X3t2e^^P‘''“i‘''“2^  )  E-^ 

(14) 

w 

II 

G.|i-‘ 

(Xitit2e^“i  -  X2rir2e^°‘M  Ei 

(15) 

^8  “  d 

(Xit3^r2e^^“i‘‘’^^^^  +  X2r  it2e^^“2'*'’^/2^ 

+  X3^X2X3e^^P‘*’°‘i'*'®2)  )  E-^ 

(16) 

We  can  use  equation  (16)  to  describe  the  overall  charac- 

teristics  of  the  device.  We  will  define  the  transmit- 

tance  of 

the  device  by; 

•'^1 


■  •'■ ’kV’-i“,  i’ 
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T  = 


and  the  net  optical  phase  shift  by; 


(17) 


if  =  tan”^ 


Im(Eg/Ei) 
Re (Eg/El) 


(18) 


IV.  SUPPLEMENTAL  PATH  DESIGN 

Two  simple  methods  for  implementing  the  sup¬ 
plemental  paths  are  shown  in  Figures  2  and  3.  The 
method  of  Figure  2  could  be  used  when  a  more  sensitive 
fiber  is  required  in  order  to  detect  changes  in  Z  or 
The  transmission  coefficient  and  phase  for  this  method  is 
described  by  the  following  equation: 


t-Ze^“z 


l-tZe^o^s 


(19) 


For  the  case  of  very  strong  coupling,  t  *  0,  the  trans¬ 
mission  coefficient  and  phase  are  completely  determined 
by  the  loop  characteristics; 


-Zeiot-  =  zei^otz+Tf) 


(20) 
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The  method  of  Figure  3  is  simply  a  way  to 
introduce  a  phase  shift  without  splicing  or  reflection 
loss.  The  transmission  coefficient  and  phase  for  this 
method  is  described  by  the  following  equation: 

E4  ir  +  Ze^“z 

—  =  - T - la—  (21) 

1-irZe^  z 

For  the  case  of  zero  coupling,  r=0,  once  again  the  output 
is  completely  determined  by  the  loop  characteristics: 


V.  DEVICE  OUTPUT  CHARACTERISTICS 

The  characteristics  described  in  this  section 
are  based  on  the  ideal  case  of  zero  coupler  loss.  Coup¬ 
ler  loss  effects  will  be  discussed  in  the  implementation 
section.  We  also  assume  that  the  device  is  under  cohe¬ 
rent  excitation  and  the  coherence  length  is  longer  than 
the  resonant  cavity  length. 

The  transmittance  of  this  device  exhibits  two 
fundamental  output  characteristics.  There  may  be  a  dip 
in  a  relatively  high  output  or  a  peak  in  a  relatively  low 
output.  These  both  occur  as  a  function  of  the  coupling 
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valuesr  transmission  coefficients  and  X2»  and  the  loop 
phase  P. 

A.  Dip  in  the  Transmittance 

This  condition  exists  when  the  device  is 
configured  to  give  a  high  off-resonance  transm ittancer  as 
shown  in  Figure  4.  When  the  proper  phase  in  the  loop  is 
achieved,  the  loop  power  builds  up  while  the  output  power 
decreases.  This  dip  is  produced  under  two  different 
device  configurations  yet  their  characteristics  are  very 
similar . 

In  the  first  case,  the  second  coupler,  CP2, 
acts  to  enhance  the  less  than  complete  coupling  of  CPl. 
For  this  condition  we  set  Xj  =  X2  =  1  and  “1  =  *  0* 

Physically,  this  represents  the  removal  of  the  supplemen¬ 
tal  paths.  Examination  of  equations  (9)  and  (16)  show 
that  resonance  occurs  for  P  =  3  '^/2,  and  under  this 
condition  the  transmittance  becomes: 

i(tir2  +  rit2  -  X3)  ^ 

T  =  -  (23) 

1  +  iX3(rit2  +  r2ti) 

At  this  point  it  is  interesting  to  relate  equation  (23) 
to  the  transmittance  of  a  single  coupler  resonator.  We 
can  do  this  by  letting  r2  *  0  and  t2  =  1  19]: 
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T 


s 


i{r2^  -  X3) 
1  +  iX3r]^ 


2 


(24) 


For  this  case  the  minimum  occurs  when  r^^  =  X3.  However, 
when  the  coupling  value  of  r^  is  not  comparable  to  X3» 
the  resonator  finesse  and  dip  magnitude  are  degraded. 
Fortunately,  the  addition  of  a  second  coupler,  as  illus¬ 
trated  in  Figure  1,  can  be  used  to  enhance  the  first 
coupler.  Referring  to  equation  (23),  the  new  condition 
for  transmittance  minimum  is: 

tj^r2  +  *^1^2  ~  ^3  (25) 

Assuming  zero  coupler  loss,  we  can  use  energy  conserva¬ 
tion  to  rewrite  equation  (25): 

-s|  1  -  Z2^  -  X3  =0  (26) 

By  inserting  the  known  values  of  r^  and  X3  we  can  find 
the  required  value  for  r2.  As  an  illustration,  consider 
the  following  example: 

NOTE:  The  loop  power  ratio  is  the  ratio  of  power  in  the 
feedback  loop  to  the  input  power. 

Single  Coupler  Resonator : 

X3  *  .9999 

r  =  .995 

t  =  .0998749 

finesse  =  617.245 

transmittance  minimum  =  .923287 

loop  power  ratio  maximum  =  383.581 
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Enhanced  Coupler  Resonator ; 

^2  =  .086  (from  equation  26) 
t2  =  .996295 
finesse  =  21,000 
transmittance  minimum  =  .000217 
loop  power  ratio  maximum  =  5,000 

This  increase  in  finesse  causes  the  magnitude  of  the 
transmittance  minimum  to  become  more  sensitive  to  the 
loss  in  the  loop.  Figures  5  and  6  are  plots  of  the 
transmittance  minimum  versus  X3  for  the  single  and  en¬ 
hanced  resonator  configurations  respectively. 

The  second  case  of  dip  in  the  output  occurs 
when  the  coupling  at  CPI  equals  the  coupling  at  CP2  such 
that : 

rj^t2  =  ^2^1  *  ^  (27) 

The  associated  device  parameters  are  as  follows:  =  X2 

=  1  and  =  02*  At  resonance  the  transmittance  is 

given  by  equation  (23).  However,  in  this  case  we  note 
that  the  transmittance  minimum  occurs  for: 

A  =  X3/2  (28) 
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Knowing  that  the  coupling  values  are  equal  and  applying 
conservation  of  energy  we  can  derive  the  following  rela¬ 
tionship  between  r  and  X3: 


Equation  (29)  shows  that  for  any  fiber  loss  we  can  find 
an  appropriate  coupling  value  which  will  minimize  the 
resonant  transmittance.  An  X3  of  .9999  yields  a  coupling 
value  of  r  =  .712  while  X3  =  .99  yields  r  =  .7553. 
Figures  7  and  8  show  the  plots  of  transmittance  minimum 
versus  X3  for  r  =  .712  and  r  *  .7553  respectively.  Un¬ 
fortunately,  the  finesse  is  still  dominated  by  the  loop 
transmission  X3.  Therefore,  even  though  we  can  optimize 
r,  the  finesse  continues  to  degrade  as  X3  decreases. 
This  is  illustrated  by  the  following  example; 

X3  =  .9999 

r  =  .712 
t  =  .7022 
finesse  =  21,000 
transmittance  minimum  =  .000328 
loop  power  ratio  maximum  =  4,998 
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r  =  .7553 
t  =  .65538 
finesse  =  165 

transmittance  minimum  =  6.343  X  10“^ 
loop  power  ratio  maximum  =  50.25 

Up  to  this  point  we  have  assumed  that  ot  = 
0^;  however#  the  transmittance  varies  when  this  is  not 
true.  For  ^  “2  location  of  the  transmittance 

minimum  changes  with  respect  to  P  as  a  function  of  1  otj^  - 
“2  !•  addition#  as  the  location  of  the  minimum 

changes#  the  magnitude  of  the  transmittance  minimum  as 
well  as  the  finesse  decreases.  This  is  illustrated  in 
Figure  9  where  is  held  constant  and  ®  2  varies  as 

shown . 


B.  Peak  in  the  Transmittance 


This  unique  condition  exists  when  the  device 
is  configured  to  give  a  low  off-resonance  transmittance. 
Figure  10  shows  a  plot  of  this  output.  The  device  para¬ 
meters  which  cause  this  output  are  listed  below: 

—  "strong"  coupling  at  CPI#  rj^  >  .9 

—  "weak"  coupling  at  CP2#  r2  <  •! 

—  Xi  *  1,  X2  »  0 


_ arbitrary 
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For  the  parameters  listed  above  and  at  resonance,  the 


transmittance  is  given  by: 


1  -  X3rit2 


Solving  this  equation  for  maximum  transmittance,  T  =  1, 


and  applying  conservation  of  energy  yields  the  following 


relationship: 


^2^  ^  "  ^1^  +  X3ri4  1  -  r22  -  1  =  0 


Unfortunately,  this  equation  only  has  real  solutions  when 


X3  «  1.  The  effect  of  X3  on  equation  (31)  is  illustrated 


by  the  following  calculations: 


X3  -  1,  ri  -  .995 


From  equation  (31),  12  »  .099878 


r2~N  1  -  rj^^  +  ^S^^l  N  ^  "  ^2^  -  1  »  0 


X3  *  .9999,  ri  *  .995 


r2  *  .099878 


•2^  1 


+  X3ri'J  1  -  r2^  -  1  *  -1  X  10"^ 
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The  previous  calculations  confirm  that  although  we  can 
minimize  equation  (31)  for  X3  l,  real  solutions  do  not 
exist  unless  X3  =  !• 

Finally,  variation  of  the  phase  just 
changes  the  location  of  the  peak  with  respect  to  P.  The 
magnitude  and  finesse  are  unaffected  by  this  variation. 

VI.  EFFECTS  OF  MODULATING  P 

Modulation  of  the  phase  shift  P  in  the  feed¬ 
back  loop  represents  a  simple  method  of  exploiting  the 
characteristics  described  in  Section  V.  Simple  sinusoi¬ 
dal  modulation,  for  example,  with  a  PZT  crystal  results 
in  a  multiple  dip/peak  output  in  the  transmittance  wave 
form.  The  general  form  of  this  modulation  is: 

P  =  Asin  (wt  +  <|)  )  +  B  (32) 

The  amplitude  coefficient  A  determines  the 
number  of  dips/peaks  in  the  output  waveform.  It  also  has 
a  small  effect  on  the  dip/peak  width  due  to  its  effect  on 
the  slope  of  the  sinusoid.  (|)  is  a  phase  bias  which 
determines  the  zero  crossing  point  of  the  sinusoid.  This 
affects  the  dip/peak  width  and  the  distance  between 
dips/peaks.  The  phase  bias  B,  which  is  external  to  the 
sinusoid,  effectively  centers  the  modulation  about  some 
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desired  point.  Changes  in  B  have  approximately  the  same 
effect  as  changes  in  <}>  • 

It  is  apparent  that  even  this  simple  modula¬ 
tion  scheme  can  have  complicated  effects  on  the  transmit¬ 
tance  waveform.  For  our  calculations  we  will  let  w  =  1, 
f  =  1/2  ffHz,  and  therefore  our  time  scale  will  be  in 
seconds.  In  the  applications  section,  the  applied  modu¬ 
lation  will  be: 

P  =  3.5  cos{t)  -  Tr/2  (33) 

This  form  is  convenient  in  that  it  is  centered  at  the 
resonance  point  and  produces  a  transmittance  waveform 
with  a  period  equal  to  seconds.  In  addition,  the 
dip/peak  resolution  is  identical  to  that  produced  under  a 
ramp  phase  modulation. 

The  Fourier  analysis  in  the  applications 
section  is  based  on  this  modulation  scheme.  A  method  for 
determining  the  spectrum  for  general  modulation  is  given 
in  the  Appendix. 

VII.  IMPLEMENTATION  CONSIDERATIONS 

This  section  contains  information  necessary 
for  the  realistic  computer  simulation  of  this  device  and 
its  actual  implementation.  In  previous  sections  two 


primary  sources  of  loss  had  been  ignored#  coupler  loss 
and  fiber  attenuation.  Data  gathered  from  the  literature 
will  be  used  to  incorporate  these  losses  into  our  calcu¬ 
lations.  The  resulting  information  will  then  be  repre¬ 
sentative  of  the  current  state  of  the  art  and  model  a 
realistic  device. 

Recent  advances  in  optical  fiber  fabrication 
technology  has  resulted  in  a  new  generation  of  low  loss 
fiber  [4,  5].  To  avoid  placing  too  severe  a  restriction 
on  the  device  design,  we  chose  a  fairly  common  fiber  loss 
value  of  .2  db/km  [13].  Choosing  a  fiber  loop  length  of 
one  meter  makes  the  loop  transmission  X3  equal  .99997697. 
Therefore  in  our  calculations  we  will  use  the  value: 

X3  =  .9999  (34) 

Coupler  loss  as  it  applied  to  a  one  coupler 
resonator  was  discussed  in  Section  II.  Since  our  device 
consists  of  two  couplers  with  different  coupling  values, 
we  will  refer  to  reference  [7]  for  the  coupler  loss 
values.  Given  the  desired  r,  the  corresponding  y  can  be 
inserted  into  equation  (1)  to  determine  the  required  t 
value.  These  numbers  will  then  be  used  in  the  transmit¬ 


tance  calculations 


When  coupler  loss  values  are  included  in  the 
energy  conservation  equation,  then  real  solutions  do  not 


exist  for  equations  (26),  (29),  and  (31).  However,  opti¬ 

mum  values  can  be  found  and  these  will  be  used  in  our 
calculations.  From  this  point  on  the  various  resonator 
configurations,  optimized  with  losses  included,  will  be 
referred  to  by  the  following  designations: 

Configuration  A  (Single  Coupler,  Transmittance  Dip) 
r  =  .995 
t  =  .09206 
Y  =  .15% 


Conf iquration  B  (Transmittance  Dip) 


n  = 

.995 

r2  * 

.09 

tl  = 

,09206 

t2  = 

.99368 

II 

1— 1 

.15% 

^2  = 

.45% 

Configuration  C 

(Transmittance 

Dip) 

^1  = 

12  =  .712 

tl  = 

t2  =  .7 

*-• 

II 

>2  =  .3% 

Configuration  D 

(Transmittance 

Peak) 

= 

.998 

r2  = 

.0632 

ti  = 

.04996 

= 

.99549 

In  Section  V  we  examined  the  characteristics 


for  configurations  A,  B,  and  C  without  losses.  The 


characteristics  with  losses 

included 

are  listed 

below: 

A 

B 

C 

finesse 

617.00 

1023.0 

976.00 

transmittance  minimum 

0.9733 

0.8754 

0.8806 

loop  power  ratio  maximum 

325.00 

0.3861 

27.000 

It  is  interesting  to  note  that  configurations  B  and  C  are 
greatly  affected  by  the  addition  of  coupler  loss.  This 
is  because  the  lower  coupling  values  required  have  a 
correspondingly  higher  coupler  loss.  Nevertheless,  they 
represent  a  significant  improvement  over  the  single  coup¬ 
ler  design.  Plots  of  the  transmittance  minimum  versus  X3 
for  configurations  A,  B,  and  C  are  shown  in  Figures  11, 
12,  and  13  respectively. 

The  most  surprising  result  of  these  calcula¬ 
tions  was  the  decline  in  loop  power  ratio  maximum  for 
configurations  B  and  C.  Although  the  finesse  is  approxi¬ 
mately  166  percent  of  configuration  A,  the  loop  power 
ratio  is  less  than  1  percent  of  configuration  A. 

Configuration  D  was  discussed  in  Section  V 
and  it  was  determined  that  even  with  zero  coupler  loss, 
real  solutions  for  r2  exist  only  when  X3  =  1*  Therefore, 
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we  assumed  a  one  centimeter  loop  length  which  yields  X3  = 
.9999995.  Comparison  of  the  idealized  calculations  to 
the  realistic  ones  are  shown  below: 


Ideal  Case  (X^  s  i) 

=  .998 
ti  =  .0632 
Yl  =  0 

finesse  =  791 

transmittance  maximum  »  1.0 
loop  power  ratio  maximum  =  250 


t2  =  .0632 


=  .998 


=  0 


Realistic  Case  (X3  =  .9999995) 
Configuration  D 
finesse  =  482 

transmittance  maximum  =  .236 
loop  power  ratio  maximum  =  59 


Therefore,  it  does  not  appear  that  this  configuration  has 
any  immediately  useful  characteristics. 


VIII.  APPLICATIONS 


It  is  obvious  that  this  device  has  a  number 
of  potential  applications.  In  this  section  we  will  exa¬ 
mine  the  use  of  configurations  B  and  C  as  microbend 
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sensors.  Potential  application  as  a  laser  gyroscope  will 


also  be  introduced. 


Microbend  sensors  are  based  on  producing  a 


change  in  transmission  along  an  optical  path.  This 


change  in  transmission  is  proportional  to  some  applied 


force  and  is  due  to  the  microbends  induced  in  the  fiber 


by  that  force.  In  reference  116],  an  analysis  of  the 


change  in  transmittance  of  a  fiber  due  to  a  single  micro¬ 


bend  is  reported.  The  study  reports  an  approximately 


quadratic  relationship  between  insertion  loss  in  dB  and 


microbend  amplitude  in  micrometers.  Applying  a  quadratic 


best  fit  to  the  tabular  data  we  can  approximate  the 


insertion  loss  for  small  distortion  amplitudes.  This 


relationship  is  given  by: 


IL  =  (2.5333  X  10"^)  DA^  +  (6.6667  x  10"^)  DA  (35) 


where:  IL  =  insertion  loss  in  dB 

DA  =  distortion  amplitude  in  urn 


If  we  assume  that  this  insertion  loss  takes  place  in  our 


feedback  loop,  then  we  can  relate  IL  to  X3  by  the  follow¬ 


ing  equation: 


IL  =  -10  log 


.9998 
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Finally,  we  can  use  equations  (36)  and  (35)  to  write  a 
quadratic  equation  in  terms  of  and  DA  only; 


Da2  +  (263.1596  x  10“^) 


DA  + 


10  log 


.9998 


(2.5333  X  10"'*) 


(37) 


Using  equation  (37),  we  can  determine  the 
relationship  between  the  distortion  amplitude  of  a  single 
microbend  in  our  feedback  loop  and  the  change  in  output 
transmittance.  Referring  to  Figures  5  and  6,  we  see  that 

I 

both  configurations  B  and  C  are  approximately  linear  over 
the  range  X3  =  .9999  to  X3  *  .999.  Tables  1  and  2  show 
the  variation  in  transmittance  versus  the  loop  transmis¬ 
sion  X3  for  both  configurations.  Using  equation  (37), 
this  range  compares  to  a  distortion  amplitude  range  of  0 
to  9.85  micrometers.  If  we  assume  a  relatively  low  input 
power  of  ImW,  then  the  output  power  ranges  are  listed 
below: 


Conf iguration  B:  .875416  mW  to  .250749  mW 
Configuration  C;  .880633  mW  to  .267967  mW 

The  sensitivity  of  this  device  is  limited 


theoretically  by  the  quantum  noise  limit  and  realisti- 
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TABLE  1.  CONFIGURATION  B 


LOOP  X3 

TRANSMITTANCE 

LOOP 

TRANSMITTANCE 

0.999900 

0.875415917980 

0.999450 

0.477098139705 

0.999882 

0.854678049932 

0.999414 

0.465406410787 

0.99986i» 

0.834421670525 

0.999414 

0.453974534365 

0.999828 

0.795304326216 

0.999378 

0.431866486022 

0.999810 

0.776419753966 

0.999360 

0.421178784610 

0.999792 

0.757969430036 

0.999342 

0.410727866912 

0.999774 

0.739942412349 

0.999324 

0.400508343438 

0.999756 

0.722328082158 

0.999306 

0.390514966571 

0.999738 

0.705116132933 

0.999288 

0.380742626246 

0.999720 

0.688296559685 

0.999270 

0.371186345779 

0.999702 

0.671859648710 

0.999252 

0.361841277835 

0.999684 

0.655795967728 

0.999234 

0.352702700547 

0.999666 

0.640096356409 

0.999216 

0.343766013752 

0.999648 

0.624751917261 

0.999198 

0.335026735368 

0.999648 

0.624751917261 

0.999198 

0.335026735368 

0.999630 

0.609754006868 

0.999180 

0.326480497894 

0.999612 

0.595094227462 

0.999162 

0.318123045017 

0.999594 

0.580764418814 

0.999144 

0.309950228349 

0.999576 

0.566756650432 

0.999126 

0.301958004261 

0.999558 

0.553063214055 

0.999108 

0.294142430830 

0.999540 

0.539676616423 

0.999090 

0.286499664887 

0.999522 

0.526589572314 

0.999072 

0.279025959156 

0.999504 

0.513794997849 

0.999054 

0.271717659498 

0.999486 

0.501286004032 

0.999036 

0.264571202237 

0.999468 

0.489055890535 

0.999018 

0.257583111578 

0.999450 

0.477098139705 

0.999000 

0.250749997107 

l.ll 
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TABLE  2 

.  CONFIGURATION 

C 

LOOP  X3 

TRANSMITTANCE 

LOOP  X3 

TRANSMITTANCE 

^ } 

0.999900 

0.880632937486 

0.999450 

0.493466010325 

0.999882 

0.960694080204 

0.999432 

0.481942031258 

»^i 

0.99986H 

0.841198048476 

0.999414 

0.470663204232 

0.999846 

0.822133964038 

0.999396 

0.459624086000 

0.999828 

0.803491262739 

0.999378 

0.448819373570 

0.999810 

0.785259683985 

0.999360 

0.438243900023 

0.999792 

0.767429260581 

0.999342 

0.427892630459 

0.999774 

0.749990308971 

0.999324 

0.417760658096 

0.999756 

0.732933419839 

0.999306 

0.407843200491 

'i' 

0.999738 

0.716249449080 

0.999288 

0.398135595891 

0.999720 

0.699929509099 

0.999270 

0.388633299708 

0.999702 

0.683964960446 

0.999252 

0.379331881108 

k 

0.999684 

0.668347403754 

0.999234 

0.370227019718 

5‘i^ 

•»y 

0.999666 

0.653068671988 

0.999216 

0.361314502439 

0.999648 

0.638120822970 

0.999198 

0.352590220365 

0.999630 

0.623496132188 

0.999180 

0.344050165800 

.5 

0.999612 

0.609187085858 

0.999162 

0.335690429379 

0.999594 

0.595186374251 

0.999144 

0.327507197277 

f. 

0.999576 

0.581486885252 

0.999126 

0.319496748510 

St 

0.999558 

0.568081698158 

0.999108 

0.311655452322 

0.999540 

0.554964077691 

0.999090 

0.303979765661 

0.999522 

0.542127468239 

0.999072 

0.296466230728 

0.999504 

0.529565488282 

0.999054 

0.289111472610 

0.999486 

0.517271925038 

0.999036 

0.281912196985 

k’K 

0.999468 

0.505240729275 

0.999018 

0.274865187902 

C»J 

ii 

1 

s 

4-’ 

I 

0.999450 

0.493466010325 

0.999000 

0.267967305627 

*t  * 

f.* 

»i^ 

■^41*, 

3 


almost  impossible  to  predict,  we  will  use  the  literature 
to  determine  a  baseline  value.  In  reference  [21]  a 
minimum  detectable  optical  power  of  3.4  x  10"^^W  was 
reported.  Using  this  as  a  guideline,  we  can  approximate 
the  sensitivity  of  our  devices.  Tables  3  and  4  list  the 
variation  in  transmittance,  for  configurations  B  and  C 
respectively,  over  the  range  X3  =  .9999  to  X3  » 

.99989999  in  increments  of  1  x  10"®.  Using  equation 
(37)  we  find  that  one  increment  corresponds  to  a  distor¬ 
tion  amplitude  of  .13  nanometers.  Assuming  the  same 
input  power  of  ImW,  the  corresponding  change  in  output 
power  is; 

Configuration  B;  P  -  1.16  nW  per  .13  nanometers 
Configuration  C;  P  =  1.12  nW  per  .13  nanometers 

These  values  are  well  within  current  detectability  ranges 
even  with  moderate  noise  in  the  system. 

We  will  not  attempt  to  relate  the  distortion 
amplitude  to  the  required  applied  force  since  that  is  a 
function  of  the  transducer  design.  However,  we  will 


point  out  that  these  values  are  for  a  single  microbend 
and  therefore  a  phenomenal  increase  in  sensitivity  can  be 
expected  when  a  typical  microbend  grating  is  used  [2,  3]. 
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TABLE  3. 

CONFIGURATION 

B 

LOOP 

TRANSMITTANCE 

LOOP 

TRANSMITTANCE 

0.999000000 

0.875675625265 

0.999899995 

0.875669832188 

0.999899999 

0.875674325733 

0.999899994 

0.875668532665 

0.999899998 

0.875673202341 

0.999899993 

0.875667409287 

0.999899997 

0.875672078955 

0.999899992 

0.875666285907 

0.999899996 

0.875670955573 

0.999899991 

0.875665162536 

0.999899995 

0.875669832188 

0.999899990 

0.875663863018 

TABLE  4. 

CONFIGURATION 

C 

LOOP  2Li 

TRANSMITTANCE 

LOOP  X^ 

TRANSMITTANCE 

0.999900000 

0.880726767376 

0.999899995 

0.880721193499 

0.999899999 

0.880725685028 

0.999899994 

0.880720111157 

0.999899998 

0.880724440529 

0.999899993 

0.880718866669 

0.999899997 

0.880723358182 

0.999899992 

0.880717784327 

0.999899996 

0.880722275838 

0.999899991 

0.880716701995 

0.999899995  0.880721193499 


0.999899990  0.880715619661 


39 


As  an  example  of  how  our  system  can  increase 
the  sensitivity  of  current  microbend  detection  sensor 
systems  we  examine  the  data  in  reference  [17].  This 
article  reports  the  design  of  a  wide-range  strain  sensor 
which  uses  a  roller  chain  to  impose  microbends  rather 
than  the  typical  corrugated  plate.  The  roller  chain  was 
used  to  increase  the  detection  range,  however  it  also 
reduced  sensitivity.  The  authors  reported  a  linear  rela¬ 
tionship  between  sensor  displacement  and  light  transmis¬ 
sion.  This  relationship  is  given  by  the  following  equa¬ 
tion: 

AT 

D  - -  (38) 

1.25 

where:  D  «  linear  displacement  in  inches 

AT  s  change  in  light  transmission. 

If  this  transducer  were  placed  in  our  one 
meter  feedback  loop,  then  equation  (38)  would  become: 

(.9999  -  X3) 

D  - -  (39) 

1.25 

Recalling  that  our  linear  range  is  approximately  X3  = 
.9999  to  X3  =  .999,  equation  (39)  gives  a  displacement 
range  of  0  to  18.3  micrometers.  If  we  assume  the  same 
input  power,  ImW,  then  the  output  power  range  is  the  same 


as  the  previous  application.  The  sensitivity  can  be 


approximated  by  referring  to  Tables  3  and  4,  and  plugging 


the  increment  value,  X3  *  1  x  10“^,  into  equation  (39). 


This  gives  a  displacement  sensitivity  of  approximately 


20.32  picometers.  Therefore,  by  placing  this  low  sensi¬ 


tivity  transducer  in  our  feedback  loop,  we  now  have  an 


ultrasensitive  linear  displacement  detector. 


Because  of  fluctuations  in  the  input  power 


level,  the  use  of  intensity  as  the  measurement  parameter 


requires  an  additional  isolated  reference  path  for  compa¬ 


rison.  Should  noise  levels  in  the  system  become  substan¬ 


tial  enough  to  affect  the  required  sensitivity,  it  may  be 


necessary  to  use  the  spectrum  of  the  transmittance  output 


for  detection.  The  magnitudes  of  the  first  10  harmonics 


for  configurations  B  and  C  are  shown  in  Figures  14  and  15 


respectively.  The  loop  modulation  is  P  =  3.5  cos(t)  - 


IT  /2.  When  we  isolate  a  single  harmonic  for  power  mea¬ 


surement  we  eliminate  a  majority  of  the  noise.  If  we 


assume  a  1  mW  power  input  then  the  power  in  the  first 


harmonic  versus  X3  is  as  shown  below; 


Configuration  B 


Configuration  C 


9999 

9998 


.01822117  mW 
.01845245  mW 


.01821135  mW 
.01843529  mW 


'S’  S'lV 
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As  a  comparison,  the  required  power  sensiti¬ 
vity  for  intensity  detection  and  harmonic  detection  for 
configurations  A,  B,  and  C  are  shown  below: 

(AX3  =  1  X  10”®,  input  power  =  1  mW) 


Configuration  A 
Configuration  B 
Configuration  C 


Intensity 
P  =  7.37  nW 

P  =  11.66  nW 
P  =  11.20  nW 


Harmonic 
P  =  17.7  pW 
P  =  25.9  pW 
P  =  25.9  pW 


Although  the  spectrum  detection  method  requires  more 
sensitivity,  it  may  be  able  to  offset  the  effects  of  a 
noisy  system. 

Another  possible  application  of  this  device 
is  as  a  laser  gyroscope.  If  we  assume  an  axis  of  rota¬ 
tion  normal  to  the  feedback  loop,  then  to  the  first  order 
the  Sagnac  phase  shift  is  given  by: 


4  ■n  NA 

'^s  - - - ^  (40) 

®  ^  C 

where:  s  “  Sagnac  phase  shift 

A  =  the  area  of  one  turn  in  the  coil 
N  =  number  of  turns  in  the  coil 
^  =  free  space  wavelength 
C  =  speed  of  light  in  free  space 
^  =  angular  velocity 
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Equation  (40)  shows  that  the  magnitude  of  the  induced 
phase  shift  is  directly  proportional  to  the  length  of  the 
loop.  Using  the  same  fiber  as  previously,  we  let  the 
length  of  the  loop  be  130.5  meters  which  corresponds  to 
an  X3  of  .997.  If  we  examine  the  transmittance  waveform 
with  P  =  3.5  cos(t)  -  11/2  over  the  range  of  t  *  0  to  t  = 
2t7,  we  notice  there  are  twc  dips;  one  is  located  at  t  = 
7r/2  and  the  other  at  t  3  ir  /2.  This  is  illustrated  in 
Figure  16.  Letting  the  loop  rotate  about  its  normal  axis 
causes  the  modulation  to  become; 

P  =  3.5  cos(t)  -  11/2  +(|)  g  (41) 

The  useful  characteristic  of  configurations  B  and  C  is 
that  the  time  between  dips  varies  as  a  function  of  ()>g. 
Listed  below  are  the  changes  in  time  between  dips  versus 


small  <l>g: 

(J)g  (rad) 

Configuration  B 

Configuration  C 

.001 

At  =  5.717  X  10"^ 

At  =  5.717  X  10"! 

.  0001 

At  =  5.65  X  10“5 

At  =  5.579  X  10"f 

.00001 

At  *  7.39  X  10"O 

At  =6.91  X  10“® 

If  we  use  a  source  with  X  =  1  micrometer  and  a  total  loop 
length  of  130.5  meters,  then  we  can  use  equation  (40)  to 
calculate  the  corresponding  rotation  rates.  If  we  assume 


3 . 5cos ( t ) 
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a  modulation  frequency  of  1  KHz,  then  we  have  the  follow¬ 
ing  relationships  between  rotation  rate  and  time  change 
between  dips: 

rad 

n  (10”®  - )  Configuration  B  Configuration  C 

s 


17.600 

At  * 

571.70 

ns 

At  = 

571.70 

ns 

1.760 

At  = 

56.50 

ns 

At  = 

55.79 

ns 

.176 

At  = 

7.39 

ns 

At  = 

6.91 

ns 

Since  inertial  rotation  sensing  requires  sen¬ 
sitivities  on  the  order  of  1  microradian  per  second,  our 
device  is  an  ideal  candidate  for  this  type  of  applica¬ 
tion.  The  major  benefit  being  the  use  of  a  single 
propagating  beam  which  makes  this  device  immune  to  the 
nonreciprocal  effects  that  plague  counterpropagating  beam 
rotation  sensors.  Additionally,  the  low  loop  power  at 
resonance  will  help  to  minimize  the  nonlinear  noise  in 
the  system  and  improve  stability. 

IX.  OPTICAL  PHASE  CHARACTERISTICS 

Up  to  this  point  we  have  dealt  only  with  the 
magnitude  characteristics  of  the  transmittance.  In  this 
section  we  briefly  describe  the  optical  phase  character¬ 
istics  that  may  bear  further  investigation. 


A1 

At  resonance,  this  device  undergoes  a 
phase  shift.  This  is  illustrated  in  figures  17  and  18 
for  configurations  B  and  C  respectively.  It  appears  that 
each  configuration  has  a  different  yet  unique  response  to 
the  phases  in  the  supplemental  paths  a  ^  and  a  2* 

Configuration  B  has  an  optical  phase  shift 
which  is  directly  proportional  to  the  phase  difference 
I  ai  -  012  I*  This  phase  shift  is  given  by  the  following 
simple  equation: 


$  *  Tr/2  +  I  ot]^  -  012  I  (42) 

If  the  supplemental  paths  are  wound  in  opposite  direc¬ 
tions  then  the  Sagnac  induced  phase  shift  in  one  path 
would  be  +<|)g  while  the  oppositely  wound  path  would  under¬ 
go  a  phase  shift  of  -  41  g.  Referring  to  equation  (42)  we 
see  that  this  configuration  could  operate  as  a  laser 
gyroscope  where  the  net  optical  phase  shift  is  given  by: 

$  *  TT/2  2  (|>g  (43) 

where:  ())g  =  Sagnac  induced  phase  shift. 

Since  the  optical  output  phase  is  directly  proportional 
to  the  rotation  induced  phase  shift,  the  sensitivity  is 
limited  only  by  our  ability  to  detect  small  changes  in 
optical  phase. 
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Configuration  C  also  has  an  optical  phase 
shift  which  is  proportional  to  the  phase  difference  |  a  ^ 
-a  2  I*  Unfortunately »  there  does  not  appear  to  be  the 
same  direct  relationship  shown  in  equation  (43).  There 
is  however  an  additional  relationship  between  the  rate  of 
change  of  the  optical  phase  shift  at  resonance  and  the 
phase  difference  ]  1  -  “2  I*  Although  the  exact  rela¬ 
tionship  is  unknown,  there  appears  to  be  a  linear  corres¬ 
pondence. 

X.  CONCLUSION 

The  two  coupler  resonator  described  in  this 
thesis  has  a  number  of  unique  characteristics.  It  pre¬ 
sents  new  ways  to  implement  old  concepts  and  introduces 
new  phenomena  with  many  potential  applications. 

State  of  the  art  single  coupler  resonators 
are  constrained  by  the  amount  of  light  that  can  be 
coupled  into  the  resonant  loop.  Our  device  presents  two 
ways  to  overcome  that  restriction  by  introducing  effec¬ 
tive  coupling  values.  These  are  listed  below  for  con¬ 
figurations  B  and  C  respectively: 

^eff  •  ^1^2  +  >^2^1  (44) 
reff  =  2rt  (45) 
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Both  configurations  outperform  the  single  coupler  reso¬ 
nator  yet  appear  to  be  equally  simple  to  fabricate. 
Should  a  coupling  technique  be  discovered  that  allows  low 
coupling  values  and  low  loss  then  the  performance  of  both 
configurations  B  and  C  will  be  phenomenal. 

As  is  true  of  all  resonators,  this  device  is 
extremely  sensitive  to  losses  in  the  loop.  That  makes  it 
particularly  useful  as  a  microbend  sensor.  In  section 
VIII  we  showed  a  typical  S/N  limited  sensitivity  to  a 
single  microbend  of  amplitude  .13  nanometers  and  in  a 
separate  case  showed  a  sensitivity  to  a  linear  displace¬ 
ment  of  20.32  picometers.  Although  these  values  were 
interpolated  from  the  literature,  they  do  indicate  the 
potential  capability  of  this  device. 

There  are  at  least  two  ways  this  device  can 
be  used  to  detect  rotation.  Both  are  immune  to  the 
nonreciprocal  effects  suffered  by  rotation  sensors  using 
counter-propagating  beams  because  they  have  low  loop 
power  maximums.  This  reduces  the  noise  effects  due  to 
Rayleigh  scattering  and  optical  Kerr  effect.  In  section 
VIII  we  discussed  a  way  to  detect  rotation  by  measuring 
the  time  between  dips.  A  rotation  rate  of  .17  micro¬ 
radians  per  second  corresponded  to  a  time  change  of  7 
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nanoseconds.  This  is  ideal  for  inertial  rotation  sen¬ 
sing.  In  section  IX  we  stated  that  the  optical  phase 
shift  in  configuration  B  can  be  made  directly  propor¬ 
tional  to  twice  the  Sagnac  phase  shift.  The  sensitivity 
of  this  method  is  limited  only  by  optical  phase  detection 


capability. 


Although  this  thesis  is  theoretical  in  na¬ 


ture,  much  of  the  modeling  is  representative  of  experi¬ 
mental  results  reported  in  the  literature.  Therefore, 
upon  fabrication,  we  should  not  expect  to  find  a  great 
deal  of  deviation  from  the  values  reported  here.  Because 
of  the  simplicity  of  the  device,  particularly  configura¬ 
tions  B  and  C,  construction  should  be  a  straight  forward 
process . 


1  A '  V,-  S’r ' 


APPENDIX 

In  general#  the  modulation  of  P  in  the  feed¬ 
back  loop  will  produce  multiple  dips/peaks  in  the  output. 
Only  in  the  special  case  used  in  this  thesis  are  there 
just  two  dips/peaks  a  distance  it  apart.  For  the  case  of 
random  distribution  the  period  is  2"^  and  the  following 
method  can  be  used  to  calculate  the  spectrum. 

1.  Calculate  the  spectrum  of  a  single  dip/peak  as  if  it 
were  alone  in  the  period. 

2.  Determine  the  time  between  the  initial  peak  and  every 
other  peak  within  the  period. 

3.  The  total  spectrum  is  calculated  by  the  following 
equation: 

Gt  “  +  e"3WT2  ^  g-jwTa  j 

where:  g(w)  »  the  spectrum  for  a  single  dip/peak. 

T  j  =  time  between  initial  dip/peak  and  the  ith 
dip/peak 
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